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We use the von Neumann entropy to study the single and many-particle on-site localizations of 
stationary states for an anisotropic Heisenberg spin-| chain. With a constructed bounded sequence 
of on-site energies for single and many-particle systems we demonstrate that the von Neumann 
entropy approaches zero indicating strong on-site localizations for all states. On the contrary, 
random on-site energy sequence does not lead to strong on-site confinement of all states. Our 
numerical results indicate that the von Neumann entropy provides a new insight to analyze on-site 
localizations for these systems. 
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I. INTRODUCTION 

Following the seminal work of P. W. Anderson in 
1958^'^, disorder induced localization became an active 
area of research in condensed matter physics.— Current 
interest in quantum computation and quantum informa- 
tion processing has further revived the importance of 
studying localization of quantum bits or qubits.^ In many 
proposed experimental realizations of a quantum com- 
puter such as ion traps^, nuclear magnetic resonance sys- 
tems^*^, optical lattices^i^, quantum dotaii, a key com- 
mon feature is the presence of qubit-qubit interaction, 
which is needed to perform two-qubit logical operations. 
The couplings of the qubits with the external environ- 
ment lead to a finite lifetime to the excited state of a 
given qubit. In order to perform a computation an ideal 
quantum computer must be isolated from the outside en- 
vironment and at the same time its qubits have to be 
accessible to perform quantum gate operations. The ma- 
jor challenge in realizing such a quantum computer is to 
obtain a delicate balance between the time for which the 
system remains quantum mechanically coherent, which 
is the decoherence time and the time it takes to perform 
quantum gate operations. 

The interactions between qubits can not be completely 
turned off in a realistic quantum computer. However, 
it can be controlled so that measurements can be per- 
formed on individual qubits. Recent theoretical works 
in this direction on many-particle systems have shown 
a way to achieve thisj ^^i^^'^** With appropriately tuning 
the on-site energy levels for the qubits it has been shown 
that delocalization can be suppressed, which results in 
strong on-site confinement of all stationary states in the 
system. In this case the effective localization length is 
much smaller than the inter-site distance instead of an 
exponential decay of the wave-function over many sites. 

Recently there has been a tremendous effort to un- 
derstand the properties of quantum entanglement, which 
is one of the most intriguing features of quantum the- 
ory. It plays a key role in many of the applications 
of quantum information processingi^ii^ Quantum entan- 
glement measured by the von Neumann entropy has 



been extensively studied in various condensed matter 
systems in recent years. It has been found that the 
von Neumann entropy can be used to identify quantum 
phase transitions in fermionic system o-'^^i^^'^^i^"'^'^^ and 
localization-delocalization transitions in interacting elec- 
tron systems22ii^. 

In this paper we have explored the possibility of using 
the von Neumann entropy as a measure for identifying 
strong on-site localization for single and many-particle 
systems of qubits. Von Neumann entropy approach pro- 
vides a new insight into localization, which is complemen- 
tary to other approaches that have been used before es- 
pecially the inverse participation ratio (IPR) li^iiii We 
have considered localization for systems with constructed 
on-site energies and random on-site energies. For a fi- 
nite system with constructed on-site energy sequence for 
qubits we show evidence for strong on-site localizations of 
all states using the von Neumann entropy. For the same 
system with random on-site energy sequence we find that 
it does not lead to strong on-site locahzation. 

II. MODEL 

Qubits being two level systems can be modeled as spin- 
i chains in a one-dimensional lattice in presence of an 
effective magnetic field along the z direction. The exci- 
tation energy of a qubit will then be the Zeeman energy 
of a spin and the qubit-qubit interaction will be repre- 
sented by the exchange spin coupling. The Hamiltonian 
of the anisotropic spin-^ chain in an effective magnetic 
field is 

H = Ho + Hi, 

L ^ L-1 
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Here Jnn+i = '^nn+i ^-iid "^nn+i the exchange cou- 
plings between nearest neighbor sites n,n + 1. L is the 
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number of sites and the parameter e„ corresponds to the 
Zeeman sphtting of spin n due to the static magnetic field 
in the z direction. We have considered Jnn+i Jfm+i 
to be equal to J so that the qubit interaction terms do not 
oscillate at qubit transition frequencies. Using Jordon- 
Wigner transformation^^ the Hamiltonian for the XX Z 
spin-chain (Eq. [T]) can be mapped onto a system of spin- 
less interacting fermions.^^ The resulting Hamiltonian, 
expressed in terms of the fermion creation and annihila- 
tion operators (a]^ and a„) is 



H — Ho + Hi, 



where 
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Hn = 



Hi 
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Here, J = Jnn+i {J > 0) is the fermion hopping inte- 
gral, JA = Jnn+i is the interaction energy of fermions 
on neighboring sites, and on-site fermion energies e„ rep- 
resent the on-site Zeeman energies of the spins (or the 
excitation energies of the qubits) . The dimensionless pa- 
rameter A (A > 0) characterizes the anisotropy of the 
spin-spin interaction. For the on-site fermion energies 
e„ we have used the bounded one-parameter sequencei^ 
with n > 1 



(-ir-^(-i; 



int(n/fc)^fc- 
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where h is the inter-subband distance that significantly 
exceeds the hopping integral J and Q!(0<a<l)isa 
one-dimensional parameter. int(n/fc) represents the inte- 
ger part. This constructed energy sequence detunes the 
on-site energies e„ from others. For a — 0, e„ splits into 
two major subbands at ±/i/2 for even and odd values of 
n. For small values of a, the two major subbands have 
width « ah and are separated by w h. With increase in a 
the subbands start overlapping and finally the separation 
between the subbands disappears. 

To compare our results with some other on-site energy 
sequence we have considered a simple random sequence 
of on-site energies 



Wr„ 



(4) 



W is the band-width and r„ are n independent random 
numbers uniformly distributed in the interval [0, 1]. 

For our model of spin-less interacting fermions there 
are two possible local states at each site, |1) and |0) cor- 
responding to the n-th site being occupied by a single 
fermion and the site being empty. The local density ma- 
trix at site n in eigenstate j is 



Pi = mH\o){o\)n + {i-m'm){i\)r. 



(5) 



The wave-function for the N-particle eigenstate is given 
by, IV'at) = J2k=i V'fcb)- In Eq-El V'™ is the amplitude of 



an energy eigenstate j at site n. 
entropy, is defined aa^ 



The local von Neumann 



Si = ~\i4f log, - (1 - log2(l - m'). (6) 

The spectrum averaged von Neumann entropy, (S) is the 
entropy averaged over all the eigenstates M and total 
number of sites N 



(S) 
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N M 
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For many-particle system with L sites and N particles, 
the number of energy eigenstates is M — L\/N\{L — N)\. 

One of the basic properties of the von Neumann en- 
tropy is sub-additivity. For two distinct quantum systems 
A and B the joint entropy for the two systems satisfy the 
inequality: S{A,B) < S{A) + S{B). The equality holds 
if and only if systems A and B are uncorrelated^, which 
is the case with different realizations of random on-site 
energies. Thus for the random on-site energy sequence 
(Eq. m we have averaged the entropies (5) over all the 
different random configurations. 

There are two special cases where Eq. [7] has simple 
analytical expressions. For an extended state j, = 
1 / Vn for all n. Then it follows that the averaged entropy 
is a constant equal to 



(S) 
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For a localized state = Snm (jm is a given site), {S) = 
0. So the normalized {S) varies between (for localized 
states) and 1 (for extended states). 

Another quantitative measure that is widely used to 
characterize localization is the inverse participation ratio, 

for eigenstate j which is defined as 



J|4 



(9) 



For our calculations we have used the inverse participa- 
tion ratio averaged over all eigenstates M 



(IPR) 



1 



M 



(10) 



For fully localized states Ij = 1. Ij « 1 is an indicator of 
strong on-site localization of the j-th state of the system. 
For extended states Ij >> 1, with no upper bound. 



III. NUMERICAL RESULTS 

A. Single Particle Localization 

We consider a single-particle system with L = 100 
and L = 300 sites. For the on-site energies we use the 
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constructed energy sequence in Eq. [3] with the following 
parameters for the inter sub-band distance h and hop- 
ping integral J : J = 1/30, ft, = 20 J. With h much 
greater than J the on-site energies are separated in two 
sub-bands for even and odd n. Next by choosing a di- 
mensionless parameter a of the order of J/h we further 
split each sub-band into two separate bands to detune the 
next nearest neighbors. We then construct the Hamilto- 
nian in Eq. [2] and diagonalize it to obtain the energy 
eigenstates. Interaction energy A is equal to zero for the 
single-particle case. With these eigenstates we numeri- 
cally calculate the spectrum averaged entropy (S) and 
the averaged inverse participation ratio (IPR) . 




Fig. [2] shows that the maximal von Neumann entropy 
'S'max and maximal IPR exhibit sharp resonant peaks 
(shown for L = 300 and L = 296 sites). Maximal en- 
tropy for each a is obtained from a set of averaged en- 
tropies for all the eigenstates, Smax = inaxj{S)j. Max- 
imal IPR is obtained similarly, IPR,„ax = maxj/j. The 
narrow peaks are seen for S'max close to zero (IPRmax 
close to one), which demonstrate that very few on-sites 
states are hybridized with each other. The resonant peak 
near a ~ 0.1 is due to hopping- induced shift of energy 
levels of « (J/2)^//i. For the chain with L — 296 sites 
the sharp peak for a « 0.1 is greatly reduced in size. In 
this case the shift in energy is « /2h. The difference 
of e„ for L — 300 and L = 296 sites is « ha^ for a < < 1 . 
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FIG. 1: Single-particle spectrum averaged von Neumann en- 
tropy (S) (normalized by the constant in Eq. |8]) is plotted for 
L = 100 and L = 300 sites as a function of a dimensionless 
parameter a. The constructed energy sequence (Eq.fSjl is used 
for the on-site energies. The hopping integral J = 1/30 and 
the inter subband distance h — 20 J is used for this plot. In 
the inset the average inverse participation ratio, (IPR) is plot- 
ted. Both (IPR) and (S) decrease with increase in a. Strong 
on-site localization (S) or (IPR) 1 occurs at a ~ 0.1 
for both L — 100 and L = 300 sites chain. 

In Fig.[T]we plot (S) as a function of a for both L — 100 
and L = 300 site chains. In the inset we show the results 
for the (IPR) . For a — > the energy sequence in Eq. [3] 
reduces to e„ — (— l)"/i/2. Hence, the stationary states 
form two bands of widths « J^/h (for h >> J) cen- 
tered at zLh/2 for even and odd n. With increase in a, 
the bands at ±h/2 further split into more and more sub- 
bands. This results in the monotonic decrease of both 
the entropy and the (IPR) . The spatial decay of single- 
particle stationary states has been studied in detail in 
Ref.^^- where it has been shown that strong on-site con- 
finement of the stationary states occur as a exceeds a cer- 
tain threshold value ath = J/2h. Our numerical results 
are in good agreement with the theoretical estimate. We 
find (5) — > and (IPR) — > 1 for a « 0.1 demonstrating 
strong on-site localization for all states. 



FIG. 2: Maximal single-particle von Neumann entropy, Smax 
(Maximal IPR in the inset) is plotted as function of a for 
L = 300 and L = 296 sites. Eq. [3] has been used for the en- 
ergy sequence. Both Smax and IPRmax sharply decrease with 
increase in a showing strong on-site single particle localiza- 
tion. The narrow peak for L = 300 sites chain close to a « 0.1 
is due to the boundary effect. This sharp peak greatly reduces 
in size for L = 296 sites chain. 

The rationale for choosing the different parameters in 
the model is shown in Fig. [3l where we have plotted the 
von Neumann entropy as a function of a and J/h for the 
single-particle system with L ~ 100 sites. The entropy 
decreases sharply with increase in a and h/ J reaching a 
minimum value (S) — 0.0031 for a — 0.166 and h/J = 
22. So for small a and h/ J ^ 20 most of the states are 
locally confined on sites. 

In Fig. [5] we have plotted the normalized probability 
distribution P{S) and P{I) (in the inset) for the ran- 
dom on-site energy sequence (Eq. U]). The probability 
distributions (histogram plots) are obtained by averag- 
ing the entropy and inverse participation ratio over 2000 
random on-site energy configurations. Sharp peak close 
to S' ^ or / ^ 1 demonstrate that many of the states 
are localized. However, the distributions are broad and 
resonances are seen indicating strong hybridization be- 
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FIG. 3: Three dimensional plot of the single particle von 
Neumann entropy, (S) as a function of a a and h/J for 
L = 100 sites and for the energy sequence. The entropy 
sharply decreases with increase in a and h/J. The minimum 
{S) = 0.0031 occurs for a = 0.166 and h/J = 22. 



tween neighboring sites, which means all the states are 
not confined. 

The entropy and the IPR for the constructed energy se- 
quence differ significantly from the random on-site energy 
sequence. For the constructed sequence we find strong 
on-site localization of all states whereas with the random 
sequence the confinement is not so strong. 



B. Many Particle Localization 

In this section we investigate the on-site localization 
of all states for a many-particle system with L — 12 
sites and = 6 particles where no two particles can 
occupy the same site. Here we have a total number of 
12!/6!6! = 924 ways of putting six particles in twelve 
sites. We construct the Hamiltonian in occupation num- 
ber basis states and then diagonalize the 924 x 924 Hamil- 
tonian matrix to obtain the energy eigenstates. 

In Fig. [5] we have plotted the entropy (and IPR in the 
inset) as a function of the dimensionless parameter a for 
the on-site energy sequence in Eq. [31 The entropy and 
the IPR have been averaged over all the 924 eigenstates. 
Three different values of the interaction parameter A = 
0, 0.5, and 1 are chosen for the plots. A = corresponds 
to only XX type coupling between the nearest neighbor 
spins - this case is similar to the single-particle case. The 
decrease in (S) with increase in a is seen for all three 
values of A. Fig. [S] shows the decrease in the maximal 
entropy and maximal IPR (in the inset) with increase in 
a for A = and 1. 

Localization of many-particle system differ from single- 
particle system due to the interaction term J A in the 
Hamiltonian. Fig. [5] shows that the entropy (and IPR) 
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FIG. 4; Normalized single-particle probability distribution of 
the spectrum averaged von Neumann entropy, P{S) is plot- 
ted for L = 300 sites with random energy sequence (Eq. |3| of 
bandwidth W — 13J. The entropy has been averaged over an 
ensemble of 2000 random on-site energy configurations. In the 
inset the probability distribution for the inverse participation 
ratio, P{I) is shown. The sharp peak at S « (and J « 1 in 
the inset) is the indication of localization. However, the dis- 
tributions are broad indicating strong hybridization between 
many sites. Multiple resonances occur due to these hybridiza- 
tion one of which we find for S ~ 0.2 or I ~ 2. 

decreases with increase in A. For A >> J/h, the energy 
bands at zLh/2 split into subbands. Such splitting lowers 
the value of (S) and IPR for non-zero A. 

For non-zero A, many particle resonances are seen in 
Fig. [5] and Fig. [6] for small values of a, indicating some 
stationary states are no longer on-site localized. Strong 
on-site localization of the states require that the energies 
of the distant sites must be tuned away from each other. 
This happens for large values of a, where we find (S) 
approaches zero as a — *■ 0.1 demonstrating strong on-site 
localization of all states. 

Similar to the single-particle case, we have calculated 
the three-dimensional plot of entropy as a function of a 
and h/J. We find that the von Neumann entropy de- 
creases with increase in a and h/J. The minimum value 
of (S) = 0.0061 occurs for a w 0.23 and h/J w 22. 

Fig.[7]shows the normalized probability distribution for 
the von Neumann entropy, P{S) (normalized probability 
distribution for the inverse participation ratio, P{I) in 
the inset) for the random on-site energy sequence. As 
before both the entropy and the IPR have been averaged 
over 2000 random on-site energy configurations. Peak at 
5 « shows localization of many of the states but not all. 
The distribution is broad and many particle resonances 
are seen (as an example the smaller peak at 5 « 0.2), 
which means that many on-site states are strongly hy- 
bridized as the stationary wave-functions spread over 
many sites. This is in contrast to the case with con- 
structed energy sequence where most of the states are 
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FIG. 5: Spectrum averaged von Neumann entropy for a sys- 
tem with A'' = 6 particles in i = 12 sites and with constructed 
energy sequence plotted against a. Three values of the inter- 
action parameter A = 0, 0.5, 1 are used for the plots. In 
the inset averaged IPR is plotted. Strong on-site localization 
(S) ^ or (IPR) 1 occurs for a close to 0.1 for all three 
values of A. For A 7^ isolated peaks at various values of 
a is due to hybridization of resonating on-site many-particle 
states. 



ied single and many-particle systems with random on-site 
energies where we find hybridization between many sites. 
All the sites for the random case are not confined which 
results in a wide probability distribution of the entropy 
with multiple resonance peaks. 

In conclusion, we have found that the von Neumann 
entropy provides a systematic approach to the prob- 
lem and is thus an effective tool to study localization- 
delocalization transitions in these many particle inter- 
acting systems. 
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strongly confined for a > 0.1. 

IV. CONCLUSIONS 

In this paper we have explored the use of von Neu- 
mann entropy to characterize strong on-site localization 
for interacting qubits. For a successful realization of 
a quantum computer strong on-site localization of all 
states is a necessity. We have considered both single and 
many-particle systems with a carefully constructed se- 
quence of on-site energies. With our numerical results 
we have demonstrated that von Neumann entropy ap- 
proaches zero for such systems with certain values of the 
parameters in the model, thus indicating strong on-site 
localization of all states. For comparison we have stud- 



FIG. 6: Plot of maximal many-particle von Neumann entropy 
(maximal IPR in the inset) for A'^ = 6 particles in L = 12 sites. 
Two values of the interaction parameter A = 0, 1 are used. 
Both Smax and IPRmax sharply decrease with increase in a 
showing strong on-site single particle localization. However, 
sharp isolated peaks are seen for A = 1.0, which indicates that 
some of the stationary states are not localized on-site. These 
states are strongly hybridized with other resonating on-site 
many-particle states. 
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